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SPECTRAL ANALYSIS OF A HIGH-VELOCITY METEOR 


Gale A. Harvey 
Langley Research Center 


SUMMARY 

A spectrogram of a fast optical meteor was reduced and analyzed, and 
60 features were identified in the spectrum. Airland ionized elements in this 
meteor radiate throughout the spectrum from 3000 A to 6800 A. A mass of 9 mg 
and an effective radiation temperature of approximately 5700 K were computed 
for the meteor. Weight ratios of Ca:Fe, Mg:Fe, and Na:Fe were computed. A 
plasma particle velocity distribution for meteors was derived, and the average 
collision speed obtained from this distribution was compared with the relative 
collision speed of a Fe-N 2 gas mixture at 5700 K. 


INTRODUCTION 

The NASA Langley Research Center (LaRC) has been active in meteor physics 
research since 1963. The research program began as a cooperative effort with 
the Smithsonian Astrophysical Observatory (SAO) and consisted of a series of 
rocket flights to simulate natural meteor conditions (ref. 1). The primary 
task of the research effort was to reduce the uncertainty in determining 
masses of optical meteors through better measurements of the meteor "luminous 
efficiency" parameter. This research effort succeeded in reducing the uncer- 
tainty in determining masses of optical meteors from 2 orders of magnitude to 
approximately a factor of 3. 

These initial studies in meteor research were not concerned with meteor- 
oid composition. To investigate meteoroid composition, the meteor spectral 
research program evolved from the NASA-SAO luminous efficiency program and used 
spectroscopy to study composition of meteoroids in the mass range from 0.001 g 
to 1 g. Most of the meteor spectra were obtained from meteor patrols in south- 
ern New Mexico. 

The NASA LaRC Faint Meteor Spectra Patrol (ref. 2) began yielding optical 
meteor spectra on a continuous basis in November 1968, and over 250 meteor 
spectra were obtained from the patrol in calendar year 1969. The patrol 
obtained spectra of faint (+i1 to -3 magnitude) optical meteors over the 
extended spectral range from 3000 A to 7000 A and had provisions for absolute 
calibration. The patrol used in-house designed and fabricated, low f-number, 
Maksutov slitless spectrographs and photoelectric meteor-detection shutter 
systems. The sensitivity of the spectrograph systems is roughly 10 times that 
of previous generation meteor spectrographs. This means that most of the 
spectra obtained by the patrol are from the free molecular flow regime and not 
from the slip, transition, or continuum flow regimes characteristic of spectra 
of similar quality obtained by earlier investigators. It should be noted that 



little quantitative spectral analysis, other than that of two very bright 
meteors (refs. 3 and ^1), had been attempted before the NASA program. 

The purpose of this paper is to report the data reduction of the first 
high-definition meteor spectrum from the NASA LaRC Faint Meteor Spectra Patrol. 
The paper then introduces a technique for the analysis of the radiative pro- 
cesses of faint optical meteors and applies it to the meteor spectrum. 

SYMBOLS 

A Einstein probability coefficient of atomic state, sec"^ 

B(T) partition function of element 

C,C constants 

E energy of atomic state above ground state, eV 

g statistical weight of atomic state 

h Planck's constant 

I line intensity 

k Boltzmann constant 

m atomic mass 

N number of atoms 

N(V) speed distribution 

T temperature , K 

« 

Tj^j meteor temperature parameter, K 

t time 

V particle speed, km/sec 

Vj, relative collision speed, km/sec 

^rms root-mean-square speed, km/sec 

Vq initial speed of meteoric particle, km/sec 

V-] lowest characteristic speed of particles for excitation, km/sec 

. o 

A wavelength, A 

o 

^max maximum wavelength, A 
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^tab wavelength of multiplet or line taken from literature 

(see table I), A 

V frequency of radiation 

- Subscripts: 

a,b designates elements (e.g., Fel) 

B Boltzmann speed distribution 

i,j ith, jth atomic states 

m meteoric particles 

A bar over a symbol denotes average. 


DATA 

A high-velocity meteor was photographed in the early morning of 
November 21, 1968, over southern New Mexico. The meteor was photographed with 
an f/0.83, 150-mm-aperture Maksutov-type slitless spectrograph. The spectro- 
graph has an inverse dispersion of 500 A/mm which is intermediate in resolution 
between the high-resolution and low-resolution spectrographs of the NASA LaRC 
Faint Meteor Spectra Patrol (ref. 2). The meteor is intermediate in brightness 
between the fainter and brighter meteors photographed by the patrol and, in 
this sense, may be considered representative of the patrol spectra. The pro- 
jected trail of the meteor passed near the Leonid radiant. The relative 
intensities from the auroral green line of 01 and from 01, NI, N2, Call, Mgll, 
Sill, and Fell confirm a velocity of 60 km/sec, or higher. 

This spectrum is quite similar to Russell’s 1958 spectrum of the upper 
trail of a -2 magnitude Perseid obtained with a 46-cm (l8-in.) Schmidt tele- 
scope with an objective prism (ref. 5). The relative strengths of multiplets 4, 
5, 20, 21, 23, and 24 of Fel; of multiplet 3 of Mgl; and of lines of Call, Mgll, 
and Sill are in general accord in the two spectra. The NASA spectrum is also 
in general accord with the lower definition faint fast spectra of Millman, Cook, 
and Hemenway (ref. 6) and Russell (ref. 7). 

The spectrum, reproduced from reference 2, is shown in figure 1. The spec- 
trum differs from previously published high-velocity spectra (refs. 8 to 10) 
of similar quality in three respects: the weakness of the H- and K-lines of 

Call, the strength and number of lines whose wavelengths are less than 3900 A, 
and the amount of air radiation in the spectrum. The H- and K-brightness 
behavior probably results from low-energy electron excitation and not from 
resonant charge exchange with N2"‘', as discussed by Hoffman (ref. 11), Harvey 
(ref. 12), and Rajchl (ref. 13). An electron density for this meteor was 
inferred from electron line densities of meteors with similar brightness and 
velocity. The lower electron density at the higher altitude of this meteor 
compared with the electron densities at the altitudes of the brighter meteors 
of references 8 to 10 probably explains the H- and K-line behavior. Merely 
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note that this behavior is not grossly anomalous in this spectrum. The 
strongest feature in the spectrum is the unresolved multiple! 3 of Mgl and 
multiple! 20 of Fel centered near 3830 A. The strength and number of lines in 
the near-ultraviolet are explained by the transmission properties of the fused 
silica of the refractive elements in the lens system, the enhanced blue sensi- 
tivity of the Kodak 2485 emulsion, and the blue blaze of the grating. Previous 
generation meteor spectrographs generally used soda- and lead-glass transmis- 
sion elements. 


ANALYSIS 

Wavelength Identifications 

The features in the spectrum are listed in table I. Sixty-three defini- 
tive features in the spectral range from 3000 A to 6800 A were measured on the 
microdensitometer tracing. Thus, the meteor spectrum is in class "a” (more 
than 49 lines) of Hillman's (ref. 14) spectral classification scheme. Identi- 
fications of the features were checked with high-dispersion spectra (refs. 10 
and 15) and with oscillator strengths (ref. 16). If a tentative identification 
was not consistent with these criteria, the feature was left unidentified. 

Three features were unidentified. The wavelengths in table I are those 

from the densitometer tracing; the wavelengths ^tab those identified, and 

are from Corliss and Bozman (ref. 16) and Moore (ref. 17). The multiplets and 
elements are from Moore (ref. 17), the products of the statistical weights and 
the Einstein probability coefficients, gA, are primarily from Corliss and 
Bozman (ref. 16), and the energy levels E-j and E2 are from Moore (ref. 17). 
Nine elements - iron, calcium, nickel, magnesium, nitrogen, oxygen, sodium, 
silicon, and hydrogen - were identified in the spectrum. 


Photometry 

The photometry is essentially that developed by Harvey (ref. I8), except 
for extension to a wider wavelength interval. The writing speed of the meteor 
image across the film was computed from the radiant position determined from 
the plate, the spectrograph elevation angle and focal length, and the meteor 
velocity at maximum light for an assumed height of 100 km. This height is 
based upon termination of the auroral green line of 01 (ref. 19) and the height 
of maximum brightness of similar meteors (ref. 20). The estimated height 
uncertainty is several kilometers and introduces an uncertainty of the order 
of 10 percent in the absolute irradiance. This is the same order as other 
errors in the photometry. 

A densitometer tracing of the meteor spectrum at the point of maximum 
brightness is shown in figure 2. The spectral irradiances of the multiplets 
of interest were obtained by reading the maximum intensity from a grid of con- 
tours of constant irradiance and multiplying by a half-width based on the plate 
resolution and the multiple! width. Thus, the irradiance values obtained are 
analogous to the equivalent widths used in studies of stellar atmospheres. The 
spectral irradiance was corrected for the emulsion sensitivity, atmospheric 
transmission, spectrograph transmission, and grating blaze efficiency. 
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Meteor Magnitude and Meteoroid Mass 

The absolute magnitude and the statistical mass of this meteor were 
obtained by comparing the spectral irradiance of the meteor in the spectral 
region from 320C A to 4800 A with the spectral irradiance of a -0.5 magnitude 
iron artificial meteor in this spectral region (ref. 21). An absolute photo- 
graphic magnitude of -1.5 (or approximately zero visual magnitude) was obtained 
for this meteor. 

Artificial meteor experiments have determined the luminous efficiency of 
iron at low meteoric velocities to within ±10 percent (ref. 1). This luminous 
efficiency may be used in conjunction with two results from the Smithsonian 
Astrophysical Observatory (SAO) - one relating the luminous efficiency of iron 
to that of natural meteors (ref. 22), and the other determining the dependence 
of the luminous efficiency upon velocity (ref. 23) - to compute the mass of 
the meteoroid. The computations yield a value of 9 mg for the mass of the 
meteoroid. 


Effective Radiation Temperature 

Spectral irradiance is proportional to line intensity; thus, photometri- 
cally corrected spectral irradiance line ratios and absolute atomic emission 
line ratios from the corresponding energy levels are equal. The line intensity 
ratio (irradiance ratio) for an optically thin gas in local thermodynamic 
equilibrium is given by 


^ g-(Ei-Ej)/kT 

Ij ’ h^jNjAj ' hVjgjAj ® 


where is the number of particles in the ith atomic state with energy Ej^ 

above the ground state, h is Planck's constant, Vj^ is the frequency of 
observed radiation from the ith state, gj^ is the statistical weight of the 
ith state, is the Einstein probability coefficient of the ith state, 

k is the Boltzmann constant, and T is the temperature. Even if the meteor 
plasma is not in local thermodynamic equilibrium, equation (1) can be used to 
compute a parameter analogous to the temperature for each choice of a pair 
of upper energy levels obtained from the observed irradiance ratios. This 
parameter can be used to study the gas. It is designated herein the "meteor 
temperature parameter T*j, and an average of T*j is the "effective meteor 
radiation temperature." 

Multiplet irradiance ratios were obtained from figure 2 and were used to 
calculate the meteor temperature parameters "^ij- Many of the spectral identi- 
fications were blends of lines; therefore, multiplets and not individual lines 
were used for the T*j determinations. Those multiplets which appeared to be 
reasonably uncontaminated by other multiplets for intensity determinations are 
listed in table II. The values for T*j calculated from the multiplets of 
table II are listed in table III. Since the magnesium lines in the meteor 
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spectrum were contaminated by iron lines, only iron multiplets were used to 
compute T*j. An effective radiation temperature of 5730 ± 670 K was obtained 
for this meteor. 

The general agreement in from different multiplets of similar 

levels, as shown in table III, indicates that assigning a temperature for 
energy levels in that range is justified. The sensitivity of the multiplet 
irradiance ratios can easily be checked via the Boltzmann formula. For an 
uncertainty in temperature of -670 K (for a temperature of 5730 K and for 
the energy levels of multiplets 5 and 41 of Fel), a corresponding mean uncer- 
tainty of 19 percent results in the irradiance ratio based on the Boltzmann 
distribution. 


Element Ratios 

The effective radiation temperature can be used to compute relative 
abundances of the dominant radiating elements with the caution that departures 
from local thermodynamic equilibrium in the plasma may introduce errors in the 
results. Other uncertainties are introduced by ionization, self-absorption, 
dissociation, and blending of lines. The Boltzmann equation can be written in 
the more general form, 


_ - -Ei/kT 

Na " Ba(T) ® 


where is the total number of atoms of the element a in a given stage of 

ionization and number of atoms in atomic state i for that stage 

of ionization; Bq(t) is the partition function of element a. The correspond- 
ing equation for line intensities is 


^a,i - h^iSa,i^a,i 




Ba(T) 


.-Ei/kT 


(3) 


for an element a. From this the number ratio of elements a and b can be 
obtained: 


^ _ ^a,i ^a^^^ Eb,jAb, j^^b, j ^_(Ej-Ei)/kT 
Nb ^b,j Eb(T) Sa,i^a,i^a,i 


( 4 ) 


By using partition functions from Corliss (ref. 24) and T^^j for T, the 
following weight ratios were calculated: 
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CaI:FeI = 0.052:1 


MgIrFel = 1.35:1 
Nal.-Fel = 0.012:1 

Because of the many uncertainties inherent in these calculations, these ratios 
are assumed to have only an order of magnitude of accuracy, and the results 
from this part of the analysis can be summarized as follows: 

CaI:FeI = 0.05:1 
MgI:FeI =1:1 
NaI:FeI = 0.01:1 


METEOR PLASMA PARTICLE SPEEDS 
Derivation of Particle Speed Distributions 

Consider a "typical" faint optical meteoroid with mass of approximately 
0.01 g and density of approximately 1 g/cm3. If this meteoroid is assumed to 
be spherical, it would have estimated linear dimensions of approximately 2.7 mm. 
At a typical meteor height of 90 km, the molecular mean free path is approxi- 
mately 2.5 cm. Thus, the free molecular flow condition is easily satisfied. 

As the meteoroid passes through the rarefied gas at this height, the air 
molecules impact the meteoroid, and as a result, either atoms or small pieces 
of meteoroid leave the parent meteoroid with virtually undiminished velocity 
relative to the atmosphere. This process continues and results in a cloud of 
meteoric atoms immediately detached from the meteoroid and moving in the air- 
stream with nearly the initial speed of the meteoroid. However, the masses of 
the individual atoms are nearly equal to the masses of the atmospheric mole- 
cules. Thus, collisions are efficient in decelerating and exciting meteoric 
atoms because of favorable momentum transfer. One may visualize a cloud of 
meteoric atoms initially moving at nearly the initial meteoroid speed. Then 
the atoms give up their energy to the atmosphere in a series of collisions 
with air molecules. For meteoroids in free molecular flow, the meteoric 
process just described is similar to the model for faint meteors formulated 
by Opik (ref. 25). 

Let Njjj(V) be the number distribution of particles of speed V in the 

1 dM 

meteor plasma, x- be the atmospheric mean free path length, and Nq = — — be 

m ds 

the number of meteoric particles entering the meteor plasma per unit distance 
along the meteor trail (M is the mass of the meteoroid). The number of mete- 
oric particles in the speed interval from Vq to Vg + Av is N^CV) Av. 
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The speed distribution can be obtained either by looking at all particles 
at one time or by following in time a set of particles as they decrease in 
speed by undergoing collisions; that is, 


Nm(V) AV = Nm(t) At (5) 

f 

where AV = Vq - V-|, At = tg - t-j, Vg is the average particle speed in the 
plasma at time tg after being sputtered by the meteoroid, and V-j is the 
average particle speed in the meteor plasma at time t-| . The number of mete- 
oric particles does not change with time; therefore, Nm(t) = Ng and 


N„(V) AV = Ng(tg - ti) 


Since tg = S./Vg and t-; = W-| for each collision, 


N„(V) AV = Ng(i 
VO 


I 


= NgS, 


-AV 


which in the limit becomes 


dV 

Nnj(V) dV = -Ng«, — (6) 

V2 


In contrast, for a gas in local thermodynamic equilibrium, the speed dis- 
tribution of particles of speed V is given by the Boltzmann distribution: 


NgCV) dV = CT-2/3v2e-C'v2/T dV 

where C and C are constants and T is the temperature. 
The average speed for the cascading system is 


Pv 


N(V) dV 


V = 




N(V) dV 


Vi 


(7) 


( 8 ) 
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Substituting equation (6) and integrating results in 


V 


m 


ViVq In (Vq/V-i) 


Vo - Vi 


( 9 ) 


where Vq is the initial meteoroid speed relative to the atmosphere and 
Vi is the lowest characteristic speed of atoms in the meteor plasma for which* 
excitation occurs. For the present example, the extreme case is assumed; that 
is, Vq = 72 km/sec (the upper limit for meteoroids) and Vi = 0.9 km/sec 
(the average speed of iron atoms at 1800 K, the melting temperature of iron). 
Then, from equation (9), 


Vjj, = 4.0 km/sec 


Similarly, the energy-averaged speed- (root-mean-square speed) for the 
cascade system can be calculated to be 

(Vrms)m = ® lo”/sec 


These speeds which may characterize high-velocity meteor plasmas can be 
compared with those for a Boltzmann distribution by using equation (7) and 
equation (8) to obtain 



Using 5700 K, the effective radiation temperature of the high-velocity 
meteor, for T and the mass of an iron atom for m gives 

Vg = 1.5 km/sec 

The relative collision speed V^, which is the_relevant speed for excitations, 
for two unlike particles with average speeds V^ and in a Boltzmann gas 

(ref. 26) is 
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For iron atoms and air molecules at 5700 K, 


Vr = ^ Vpg = 2.8 km/sec 


Thus, the average or characteristic collision speeds of particles of a 
meteor plasma derived from an inverse square velocity distribution of meteoric 
atoms (Vqj) generally agree with the average or characteristic speeds obtained 
froitf the observed spectra by assuming a Boltzmann speed distribution (Vj,-). 
These meteor plasma characteristic particle speeds are more than an order of 
magnitude lower than the initial meteoroid speed. 


Collision Speed and Radiation 

From table I in conjunction with figure 2, it can be seen that, in general, 
the low energy levels (3.0 to 4.5 eV) of the atoms are heavily populated and 
intermediate and higher levels (greater than 4.5 eV) are sparsely populated. 

For example, multiplet 5 of Fel in the wavelength range from 3720 A to 3750 A 
is very strong, but intrinsically strong lines at 3686 A, 3694 A, and 3701 A 
of multiplets 386 and 394 are not detected in the spectrum. The upper energy 
level of multiplet 5 is 3.3 eV. The upper energy level of multiplets 386 
and 394 is 6.3 eV. Most plasmas tend to emit most of their radiation at ener- 
gies almost an order of magnitude greater than the characteristic speeds of the 
particles. This generally observed phenomenon is predicted for blackbody radi- 
ation by Wein's displacement law (^max”^ = Constant) and the definitions of most 
probable, average, or rms speeds (V » 2kT/m) . The cascade model gives a rms 
speed for the cascade of atoms from meteoroid speed, 72 km/sec, to 0.9 km/sec 
of approximately 8 km/sec. The cascade velocity distribution could explain 
the general lack of extremely high (greater than 10 eV) excitation and ioniza- 
tion of meteoric atoms (i.e., excitation and ionization comparable with the 
meteoroid velocity). However, even the collision-driven speed distribution 
of the cascading system results in a higher average relative particle speed, 

4.0 km/sec, than that inferred from the observed spectra, 2.8 km/sec. Thus, 
the actual meteor spectra may indicate that the higher energy levels of the 
atoms are less populated than may be predicted by a simple velocity cascading 
system. This may result because some of the radiation may not occur during the 
main downward cascade in velocity of the meteoric atoms but rather during the 
relaxation of the heated atmospheric-meteoric plasma after passage of the mete- 
oroid. Thus, the radiation may be even more nearly in thermodynamic equilib- 
rium than indicated by cascade excitation. This is consistent with the strength 
and time behavior of the ground-state transitions observed in chopped meteor 
spectra and meteor flares (ref. 27). 


CONCLUDING REMARKS 

Analysis of the spectrogram of a high-velocity optical meteor indicates 
that although the initial meteoroid entry velocity was 60 km/sec or higher 
(allowing maximum excitation energies over 500 eV), most of the radiation from 
the meteor is due to transitions from less energetic states (3.0 to 4.5 eV) . 
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This implies that the colliding particles responsible for the excitations have 
a speed distribution characterized by average speeds considerably less than the 
initial meteoroid velocity. This in turn is an indication that the meteoric 
plasma may, for some purposes, be approaching conditions of local thermodynamic 
equilibrium and that the Boltzmann relationships may be profitably used to 
gain insight into the conditions of the meteoric plasma and the elemental abun- 
dances of the meteoroid. 

A comparison of the particle speed distribution predicted by a simple 
velocity cascade model for meteor atoms deceleration with the Boltzmann speed 
distribution provides an indication of the degree of local thermodynamic equi- 
librium for the meteoric plasma. The average relative collision speed pre- 
dicted by the model is ^.0 km/sec which compares with 2.8 km/sec for average 
relative collision speed derived from the observed spectrum with a Boltzmann 
equation. 

Using the Boltzmann relationships, an "effective radiation temperature" 
of 5730 ± 670 K was calculated for the high-velocity meteor on the basis of 
ratios of Fel multiplet transitions identified in the spectrum. 

The Boltzmann relationship was also used to calculate relative weight 
ratios of neutral elements identified in the spectrum. To an estimated order 
of magnitude accuracy, the results of the calculations are 

Cal: Fel = 0.05:1 

MgI:FeI =1:1 

NaI:FeI = 0.01:1 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23666 
May 31, 1977 
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Multiplet and 
element 


V, 

A 

^tab, 

A 

Multiplet a 
element 

3020 

3020.5 

9 

Fel 


3020,6 

9 

Fel 


3021 . 1 

9 

Fel 

3110 

3099.9 

28 

Fel 


3100.3 
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TABLE II.- PRINCIPAL MULTIPLETS IN METEOR SPECTRUM 


Element 

Multiplet 
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3.35 
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TABLE III.- T*j 


CALCULATED FROM EQUATION (1) USING 
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parameter, T*j, 
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6,23 
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Average . . . 


5727 
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Figure 1.- Spectrum of high-velocity meteor photographed on November 21, 1968, 

over southern New Mexico. 





Contour lines represent constant 
irradiance in 10"° ergs-cm^-A-sec at 
a distance of 100 km from the source 



NASA-Langley, 1977 
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WAVELENGTH, A 

Figure 2.- Densitometer tracing of meteor spectrum from position A (fig. 1) along trail. 
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